Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences

The study of tissue dynamics has been stimulated during the last decades thanks to the use of quantitative descriptions, with the development of several theoretical and computational frameworks, many of them revolving around the notion of reaction-diffusion systems, eventually with additional structure variables beyond time and space. The use of structure variables can accommodate phenotypic traits. In this work, we study a family of competition models, where a given population depends on a resource (e.g. oxygen) and several populations are competing for it. Our quantitative description incorporates phenotypic traits and heterogeneity at the level of cell cycle variations, which influence replication rates via oxygen consumption. This enables us to replicate the fitness of specific subpopulations to environmental conditions (e.g. oxygen shortage or external influences). Using numerical simulations, we show that such models display dynamical pattern formation in the form of coupled travelling wave profiles that expand or retreat at the same wave speed. The full theoretical analysis of such dynamics is quite involved; to circumvent this difficulty, we introduce a quasi-stationary approximation for the resource dynamics. We find that this approximation can reproduce the overall behaviour very accurately, with the additional benefit of allowing theoretical treatment of the reduced model. In this way, we provide estimates on the wave speed which are numerically shown to be robust across a wide range of macroscopic parameters of the full model. The wave speeds are thus found to depend strongly on the proliferation rate of the fittest population, resembling a winner-takes-all dynamics.

BackgroundPlatelet activation in flowing blood is a multiscale process in which vessel-scale hemodynamics, red blood cell (RBC) mechanics, adhesive receptor interactions, and intracellular signalling jointly determine thrombotic risk. Individual components are well studied, but a single reduced description that carries each explicitly from vessel-scale flow to mechanosensitive calcium entry, with dimensionally consistent couplings, remains uncommon. ObjectivesWe develop and analyse a reduced, six-module mechanobiological framework for platelet priming spanning the cascade from hemodynamic shear to mechanosensitive calcium entry, and we delineate which elements are supported by existing evidence and which are new, testable hypotheses. MethodsThe framework comprises six coupled modules: (I) hemodynamic forcing from the incompressible Navier-Stokes equations, with an objective principal-strain-rate measure for extensional flow; (II) RBC-mediated platelet margination and near-wall delivery, closed by a near-wall arrival flux; (III) von Willebrand factor (VWF) activation with a bounded kernel and glycoprotein Ib (GPIb) catch-slip capture, resolved through an explicit contact area and a bond-dependent mobility that progressively immobilises wall-interacting platelets; (IV) a single-load membrane-stimulus formulation; (V) mechanosensitive gating and a dimensionally consistent cytosol-store calcium model with extracellular influx; and (VI) a phenomenological mechanical-memory state. We formally derive that the single-platelet stochastic dynamics and the continuum population balance form a Fokker-Planck pair, with the spatially varying diffusivity handled by an explicit drift correction. ResultsThe framework yields a family of mechanochemical dimensionless groups delineating priming regimes. Its central prediction is reformulated as a falsifiable, history-sensitive signature: in a conditioning-test protocol, a low-tension conditioning block charges the memory state, and a fixed sub-threshold test pulse then reports a delay-dependent calcium facilitation that decays on the memory time{tau} m and is distinguishable from no-memory gating, channel adaptation, and residual-calcium priming. We show explicitly that the previously proposed pulsatile-versus-monotone contrast is a nonlinear convexity/thresholding effect of the gating nonlinearity--its difference-in-differences is approximately zero-- and is therefore not a valid test of memory; the conditioning-test signature is. A second prediction links RBC stiffening to reduced near-wall delivery and captured-platelet calcium response, upstream of intrinsic platelet signalling. ConclusionsThe framework provides a dimensionally consistent, mechanistically grounded and hypothesis-generating description linking hemodynamic forcing to mechanosensitive calcium entry. It demonstrates how history-dependent platelet priming may arise from a phenomenological sensitisation state and proposes a conditioning-test protocol for comparison against adhesive, channel and intracellular-store persistence. The framework is calibratable rather than validated, and the quantitative outputs shown use representative uncalibrated parameters.

Structural network representations of metastatic dissemination typically focus on static topology without resolving transport dynamics, relaxation timescales, or steady-state behaviour. Here, we formulate a discrete Markovian transport model on a directed higher-order network with transition rates derived from qualitative clinical affinity classes. By constructing a non-Hermitian row-stochastic transfer operator, we characterise the relaxation dynamics through its spectral decomposition. The system exhibits a fast-mixing regime characterised by a spectral gap of {gamma} {approx} 0.67, corresponding to a characteristic relaxation timescale of {tau} {approx} 1.49 discrete steps, with the influence of the primary tumour origin progressively attenuated during dissemination. Convergence towards a non-equilibrium steady state (NESS) is accompanied by a reduction in Shannon entropy, concentrating probability mass within specific topological sinks. This spectral relaxation delineates two distinct dynamical regimes: early transient dissemination (n < {tau}), dominated by local organ-specific transition probabilities (organotropism), and the asymptotic regime (n > {tau}), determined increasingly by the global transport architecture of the network. Comparison with independent clinical and autopsy observations across 21 primary tumours and 23 target organs indicates that the predicted stationary distribution is consistent with the observed hierarchy of metastatic organ involvement.

Mechanical properties of biological tissues, driven by passive and active forces, play a vital role in several processes ranging from development to cancer metastasis. However, the dynamical responses of cells in tissues, subject to mechanical deformations such as shear and the associated rheological properties, are not well characterized. Here, we use three-dimensional agent-based models for normal and cancer tissues to investigate their responses to simple shear as a function of cell stiffness and stochastic active forces. In the normal epithelium, with uniform strength of active force, the yield stress as a function of shear rate follows the Herschel-Bulkley form over a range of cell volume fraction. Strikingly, the shear rate dependence and the elasticity-dependent changes in the yield stress fall on master curves upon suitable scaling. To model cancer-like behavior, a certain fraction (Np) of cells was chosen to have enhanced activity and decreased stiffness. As Np increases, the extent of collective cell movement decreases, transitioning from affine (collective) to non-affine (individualistic) movement, a finding that is in accord with imaging experiments. Simulations of a model of a stiff solid tumor, with radius Rs embedded in normal tissue, show that as Rs increases, the yield stress increases. Interestingly, the cells migrate collectively as Rs increases. A Gaussian Mixture Model (GMM) and a mean field theory quantitatively account for the simulation as well as experimental results on cancerous, non-cancerous, and a mixture of these two types. The combined theoretical and experimental study establishes that heterogeneity in stiffness and activity determines non-affine movements in normal and cancer tissues.

A fundamental biological process of transcription occurs in the cell nucleus, which is a complex medium that also contains multiple heterogeneous structures known as biomolecular condensates. Interestingly, some of these condensates contain mRNA molecules in addition to proteins, suggesting an important cellular role in transcription that is not yet well understood. In this work, we develop a minimal theoretical framework for quantitative investigation of the role of reversible mRNA condensation in transcription. Our discrete-state stochastic approach accounts for the most relevant processes, allowing us to explicitly evaluate the properties of the system and clarify the effects of condensation. Analytical calculations supported by computer simulations suggest that reversible mRNA condensation influences the transcription processes by maintaining a constant level of free mRNA in the nucleoplasm while lowering the degree of stochastic noise and increasing the robustness against external perturbations. Physicochemical arguments are presented to explain these observations. The proposed theoretical framework elucidates important microscopic aspects of transcription, providing a convenient quantitative tool for investigating complex biological phenomena.

Tumor spheroid mechanics arise from both the mechanical properties of individual cells and the adhesive interactions that organize them into tissues. The relative contribution of these two factors to the bulk mechanical behavior, however, remains difficult to disentangle experimentally. Here, we develop a computational model of micropipette aspiration to compare the mechanical response of isolated cells and multicellular spheroids within a common computational framework. By independently varying single-cell stiffness and cell-cell adhesion, we quantify their effects on aspiration dynamics, effective elastic modulus, and viscoelastic relaxation. Our results show that increasing single-cell stiffness substantially alters the mechanics of isolated cells but has limited influence on the effective elastic modulus of multicellular spheroids. In contrast, changes in cell-cell adhesion produce pronounced effects on spheroid effective elastic modulus. Nevertheless, both parameters increase the retardation time governing the transition from the initial elastic response to long-time viscous deformation. These findings suggest that multicellular elasticity is governed primarily by intercellular mechanical coupling, whereas the dynamical response to applied stress depends jointly on cell-scale mechanics and cell-cell adhesion.

Contact inhibition of proliferation (CIP) enables epithelial tissues to self-regulate growth and maintain tissue homeostasis. However, how cell-level mechanical contact, tissue-scale structural order, and proliferation kinetics interplay remains a fundamental open question in living matter physics. Here, we present a particle-based model of a confluent epithelial monolayer governed by overdamped dynamics, where individual cells interact via a two-dimensional hard core- soft shoulder potential. By comparing structural evolution during quasistatic densification with previously reported experimental division kinetics, we find that the dynamics of proliferation arrest mimics the onset of direct steric contacts between the hard cores of the shell. Identifying hard core contacts as the physical driver of CIP, we couple our mechanical model with a stochastic Monte Carlo division scheme in which the instantaneous division rate decreases to zero from an intrinsic value as the number of hard core contact increases to six from zero. We demonstrate that for high intrinsic division rates, the cellular densification outpaces mechanical relaxation. This kinetic mismatch drives premature hard-core contact formation, shifts the onset of jamming and contact inhibition to lower packing fractions, and induces increasingly disordered transient configurations before the tissue universally converges to a hexagonal close-packed limit. Our model's predicted division kinetics and structural order evolution are consistent with epithelial monolayer experiments, both reported and our own. This minimal physical framework links single-cell steric contact mechanics directly to tissue-scale growth regulation and structural evolution.

Microbial populations generally evolve in fluctuating environments under time-varying conditions. These are often described by binary switching models, sometimes seen as coarse-grained feast-famine cycles, in which resource availability switches abruptly between abundant and scarce conditions. However, experimental studies suggest that feast-famine environments actually exhibit more complex temporal dynamics. Here, we study how two strains, one growing slightly slower than the other, compete for the same resources in fluctuating environments comprising a finite number of intermediate states, each having its own carrying capacity. Environmental switching between these states and their carrying capacities represents gradual changes in nutrient availability. This class of multi-state stochastic switching models can be interpreted as a coarse-grained description of feast-famine cycles and allows us to investigate strain competition under the gradual recovery and depletion of resources. By computational and analytical means, we characterise the population dynamics in these multi-state fluctuating environments. In particular, we study how the switching rates and distribution of carrying capacities affect the population-size statistics, fixation probability, and mean fixation time. By comparing these results with their counterparts in binary environments, we clarify how the frequency and amplitude of environmental fluctuations influence population dynamics in coarse-grained feast-famine cycles.

A classical study found no excitation transfer when isolated cardiomyocytes were placed side by side, whereas a recent paper reported action potential transfer in carrdiomyocytes placed end to end. We use nanoscale numerical simulations based on the full Poisson-Nernst-Planck equations to investigate whether these apparently opposing observations can be explained by the different geometrical configurations. The computations show that in the end-to-end configuration, ephaptic coupling occurs when the intercellular cleft is sufficiently narrow and a sufficiently large fraction of the sodium channels is localized at the intercalated disc. Coupling is strengthened when the sodium channels are concentrated in fewer clusters and when ionic diffusion within the cleft is reduced. Under these conditions, excitation transfer occurs on a timescale consistent with rapid cell-to-cell activation. Conduction depends biphasically on cleft width and terminates abruptly beyond a critical width. Localization of potassium channels at the intercalated disc has only a moderate effect, whereas gap junctions substantially improve conduction and reduce the relative contribution of ephaptic coupling. In the side-by-side configuration, excitation transfer does not occur under physiological conditions and requires highly flattened cells, minimal separation, and unrealistically strong sodium-channel clustering. The different outcomes of the side-by-side and end-to-end experiments can therefore be explained by the fundamentally different geometrical conditions for ephaptic coupling.

Within the studies of stochastic gene expression, apart from the variability of copy number of gene products, the problems of threshold crossing of those products are biologically important as they often lead to terminal cellular events. Here, we study the threshold crossing problem of the messenger ribonucleic acid (mRNA) and present an exact probability distribution of first passage times in Laplace space. The function furnishes moments of any order and also predicts the characteristic time of the exponential tail of the distribution, which we match against Gillespie simulations. We find that all the measures of relative fluctuations of the threshold crossing times show U-shapes within this simple model of mRNA, as was found earlier in more mathematically involved models of threshold crossing time statistics of proteins. Furthermore, we extend the exact formula to include the phenomenon of DNA duplication and the corresponding doubling of transcription rate. As expected, the distribution varies considerably depending on the onset of the duplication stage within the cell cycle.

Purpose Prevention and early detection of osteoporosis remains a global challenge, more so in regions like the Philippines where screening barriers exist. Chest x-rays meanwhile are relatively inexpensive, and more frequently done, and therefore can be used for opportunistic screening. This study aimed to develop a deep learning model for osteoporosis detection from chest x-rays using DXA as the gold standard. Methods A convolutional neural network called Osteo-AI was developed using 406 pairs of chest x-rays and DXA scans of Filipino patients aged 50 and above. With data augmentation, the training set expanded to 6,300 pairs. Gradient-weighted class activation mapping technique was applied to localize and identify patterns and areas in the chest x-ray images correlating with osteoporosis. Results Training data consisted of 369 female patients and 37 males. Ages of the patients ranged from 50 to 89 with a mean age of 63 years old. Initial testing yielded promising results, with Osteo-AI achieving a diagnostic accuracy of 85.71%, easily outperforming a benchmark of 33.33% Conclusion Our findings suggest the potential of Osteo-AI to enhance osteoporosis screening accessibility, aiding in early intervention to prevent fragility fractures. Further research involving larger datasets is warranted to refine and optimize the model, potentially improving detection accuracy and expanding its utility in global healthcare settings.

Maternal thermal strain is associated with adverse pregnancy outcomes, yet fetal temperatures cannot currently be directly measured, limiting quantification of fetal thermal strain. Here, we develop two steady-state models for estimating internal temperatures in a near-term fetus. First, we improve the only previously published human fetal thermoregulation model, deriving a closed-form solution within its simplified uniform-cylinder representation. Second, we introduce a multilayer, anatomically segmented model that resolves tissue-specific temperatures. Both couple the fetal body to central blood pool and amniotic fluid compartments and incorporate a new placenta-umbilical cord heat-exchanger representation. Predictions agree with available intrauterine scalp measurements, with fetal core and head-center temperatures approximately 0.5{degrees}C and 0.8{degrees}C above maternal core, respectively. Physiologically plausible changes in umbilical cord heat-exchanger effectiveness or blood flow increased fetal temperatures by approximately 0.3{degrees}C. These models enable estimation of otherwise inaccessible temperatures, while the multilayer formulation lays a foundation for transient, coupled maternal-fetal thermoregulation modeling.

Classical neurogeometric models describe the primary visual cortex as a fibered structure in which retinal position and local orientation are coupled through the geometry of the roto-translation group. We extend this approach to the visuomotor cortex by modeling it as an assemblage of visual and motor cortical geometries. The model combines orientation-selective representations, analogous to those of the primary visual cortex, with movement-direction-selective representations, analogous to those of the primary motor cortex, in order to describe the mixed visual and motor selectivity observed in the visuomotor cortex. We introduce a coupled visuomotor structure in which visual orientation and motor direction coexist over a common spatial plane and interact through a relative-orientation constraint. Neural responses are modeled by orientation- and direction-dependent profile functions, and preference maps are obtained from vectorized population responses. Numerical simulations generate visual, motor, and mixed visuomotor response maps. A competition rule between visual and motor responses produces incidence ratios close to experimental observations in macaque visuomotor cortex. This framework provides a first neurogeometric approximation of visuomotor functional architecture and a mathematical setting for studying visually guided action.

The development of neural connections in the brain results from a complex interplay between biological processes and mechanical forces. A key question in neuroscience is how physical forces and the mechanical properties of brain tissue influence the formation of structural connections. Here, we demonstrate that mechanical forces play an essential role in shaping the emergence of short-range connections, particularly U-shaped fibers that link neighboring regions of the cortex. Using a computational model that incorporates our "stress-dependent axon reorientation" hypothesis, we simulate how growing axons respond to the mechanical stress field generated by cortical folding. Our results suggest that axonal growth and reorientation may be strongly influenced by local mechanical cues, helping establish the organization of these short-range pathways. Supported by in vivo diffusion tensor imaging and histological observations, our findings provide a physical explanation for why these fibers predominantly adopt U-shaped trajectories, and why connections between gyri (ridges) are more prevalent than those between sulci (valleys) or spanning gyri and sulci. These results suggest that understanding the mechanics of brain folding is critical for fully explaining the formation of brain connectivity and its variations in health and disorder. Teaser: Mechanical forces during cortical folding guide the formation of short association fibers in the brain.

This study documents the derivation and validation of a deterministic algorithm for cause-of-death (COD) ascertainment from longitudinal real-world medical claims data, evaluated against an independent state-level death certificate file. Death certificates are the dominant reference standard in mortality research but carry well-documented limitations, including primary-cause error rates estimated at 20-40\% across empirical studies. A matched analytic cohort of 216,382 individuals (Connecticut death records, 2017--2025, age 25 and above) was constructed after exclusion of mechanism-of-injury cases and removal of ill-defined symptom-code entries from both sources. Concordance between algorithmic and certificate-based COD was assessed through three complementary frameworks: age-stratified positive predictive value (PPV) at the ICD-10-CM chapter level under a full-set concordance scenario; mean absolute rank difference (MARD) for chapters identified by both sources; and analyses of breadth, depth, and code-level specificity of COD reporting. Chapter-level PPV was strongest for individuals aged 55 and above, with all estimates representing conservative lower bounds given the known error rate of the certificate reference standard. The algorithm consistently reported broader and more granular contributing cause profiles than the death certificate, with discordances directionally consistent with the well-documented tendency of certificates to under-report contributing conditions. These findings support the conclusion that algorithmic COD ascertainment from longitudinal claims data is a feasible and scalable alternative to certificate-based attribution and, at population scale, a principled methodology for characterising death certificate error rates beyond what small-sample chart review studies can achieve.

Wound repair is an important mechanism to preserve tissue integrity in organisms after injury. However, why different tissues exhibit different mechanisms to repair wounds is a long-standing question that remains unanswered. In this work, we theoretically explore the role of the purse string, an actomyosin contractile cable used by tissues to close small wounds. Does the tissue 3D geometry influence the efficiency of the purse string in driving wound closure? Using a 3D biophysical model, we study in silico tissues with the same cell volumes but different aspect ratios, ranging from squamous to thick and tall tissues. The model predicts that taller cells are easily deformed by the purse string. In contrast, very squamous cells require a very strong purse string that might demand additional cellular mechanisms to close the gap. These findings establish a theoretical framework to predict the optimal biophysical mechanisms of wound healing in different tissues. Graphical abstractCells of different aspect ratios can be observed in a range of organisms with different function and mechanics. The wound healing efficiency of the purse string increases with the cell aspect ratio in our theoretical exploration. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=135 SRC="FIGDIR/small/743165v1_ufig1.gif" ALT="Figure 1"> View larger version (23K): org.highwire.dtl.DTLVardef@d44ab0org.highwire.dtl.DTLVardef@1737cbaorg.highwire.dtl.DTLVardef@101b5d4org.highwire.dtl.DTLVardef@1487f26_HPS_FORMAT_FIGEXP M_FIG C_FIG

Within a single human lifetime, the spectral environment has been fundamentally reshaped. Broadband daylight, rich in infrared (IR) photons arising from solar and atmospheric physics, has been replaced in the built environment by narrow, engineered spectra that largely exclude long wavelengths. While modern lighting is optimised for vision, the non-visual photobiology of metabolism may depend on spectral components that are now absent. When expressed in photon-energy units, the solar spectrum exhibits a broad maximum near 0.75 eV. This range overlaps with the activation and reorganisation energies governing mitochondrial electron-transfer kinetics. Within a Marcus-type framework, IR photons are therefore positioned to modulate rate-limiting metabolic steps by biasing barrier-crossing probabilities rather than supplying chemical energy. These wavelengths also penetrate deeply into tissue in a scattering-dominated regime, forming a diffuse internal photon field capable of interacting with distributed mitochondrial networks. We propose the term photometabolism: a solar-driven, non-photosynthetic modulation of core metabolic processes. A scaling analysis shows that photon interception in this band varies with body mass in parallel with basal metabolic rate, suggesting that ambient sunlight provides sufficient flux to influence metabolic kinetics across the biosphere. These findings have implications for physiology, ecology and the design of indoor environments whose lighting spectra increasingly diverge from their evolutionary context.

The collective spatial phenomenon of complement protein opsonization on nanoparticle surfaces is a key component of the immune response to viruses, engineered nanoparticles, and diseased cells. Recent work showed this opsonization follows a sharp, percolation-like transition versus the spacing d between surface-bound attachment sites, leaving two open questions: 1) whether the transition exhibits hallmarks of true criticality, such as diverging susceptibility, and 2) whether it can be distinguished from an alternative first-order cooperative (Hill-type) process producing an equally sharp threshold without true criticality. Here, we resolve both questions using a hierarchical statistical-mechanics treatment spanning stochastic, mean-field, and spatial reaction-diffusion models. The variance of two order parameters, peak complement activity and activation lifetime, diverges near threshold and sharpens systematically with system size, the defining signature of a critical point rather than a smooth cooperative response. Extending the analysis across site spacing and intrinsic kinetic rate constants traces a two-dimensional locus of critical points with consistent critical exponents throughout, establishing a single, robust universality class. The mean-field dynamic exponent for activation lifetime agrees quantitatively with the exact value predicted for the general epidemic process. Finally, a reaction-diffusion model of the nanoparticle surface shows the critical locus is set by a diffusion-limited length scale, establishing complement percolation as a fundamentally transport-limited surface reaction. These results place complement activation within the percolation universality class and identify the physical parameters, diffusion, catalysis, and decay, that govern its critical threshold, with direct implications for rational design of complement-evading nanomaterials, immunology, and evolutionary biology.

Time compaction theory is a general framework explaining how a brain can efficiently deal with dynamic situations occurring in, e.g., sports games. It involves a geometric representation of the time dimension, which enables effective learning and strategic action planning. The theory has recently received experimental support in humans. However, its current computational model has an important limitation: it does not account for deliberate waiting and speed modulation, behaviors ubiquitous in natural environments. This work substantially extends the original model formulation by a dimensional lifting of an n-D workspace into (n + 1)-D mental space, where time remains geometrically embedded. The proposed biologically inspired computational model can generate adaptive behavior across increasingly complex situations, from navigation in everyday social environments to competitive sports. Furthermore, by actively conditioning the expected responses of other agents and stabilizing future predictions, we introduce the concept of uncertainty points in sequences of generalized cognitive maps to support the generation of adaptive strategies in interactive environments, where future prediction has a limited time horizon. Thus, we provide a mechanism for chaining short-term solutions into long-term strategies, which is illustrated by simulating the behavior of a player in a real football game. Author summaryHumans often anticipate future interactions in dynamic environments. Many behaviors, such as avoiding other pedestrians, letting someone pass through a narrow corridor, or reproducing the kind of dribbling maneuvers performed by elite football players, require deciding not only where to move but also when to move. Existing theories suggest that the brain simplifies such situations by representing future interactions as static spatial maps, making them easier to learn and recall. However, current computational models cannot naturally account for common behaviors such as waiting, slowing down, or modulating speed. Here we show that these behaviors readily emerge if the model space is extended by an additional virtual coordinate that encodes accumulated waiting rather than physical time. The proposed model simultaneously admits a wide variety of behaviors, including speed modulation, multigoal decisions, and compound actions, while preserving the principles of time compaction. We illustrate the model in everyday situations and by reproducing two real football plays, comparing the observed behaviors with model simulations. Our results suggest computational principles through which the human brain may efficiently represent, memorize, and exploit dynamic situations.

Recent experiments show that protein condensates sitting on opposite surfaces of a flat lipid membrane move together and prefer to overlap, even though they cannot touch each other. This points to an indirect, membrane-mediated interaction. Two mechanisms could be responsible: a curvature-induced interaction, which is energetic in origin, and a fluctuation-induced interaction, which is entropic. Here we study both with coarse-grained molecular dynamics simulations, using Cookes implicit-solvent lipid model together with a generic bead-spring polymer model for the condensate. We compute the potential of mean force between two condensates across the membrane. For condensates of the same size, full overlap is unfavorable, and the pair instead settles into a partially overlapping state that bends the membrane into an S-like shape. When the two condensates differ strongly in size, full overlap becomes favorable. We explain this with a simple geometric picture. The condensate wets the membrane as a thin film and imposes curvature only along its rim, while membrane tension flattens the membrane under its interior. The resulting ring of curvature can trap a smaller condensate on the opposite side. We also compare the bending undulations and the effective bending modulus of a bare membrane, a membrane with one condensate, and a membrane with condensates on both sides. A wetting condensate suppresses the undulation modes and stiffens the membrane, but whether this makes overlap entropically favorable remains inconclusive. Our results indicate that the coupling is driven mainly by curvature, and that it depends on the wetting mechanism and on the membrane tension.